Perfect State Transfer in T'T-symmetric Non-Hermitian Networks 
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We systematically study the VT (parity-time reversal) -symmetric non-Hermitian version of a 
quantum network proposed in the work of Christandl et al. [Phys. Rev. Lett. 92, 187902 (2004)]. 
The exclusive nature of this model show that it is a nice paradigm to demonstrate the complex 
quantum mechanics theory for the relationship between the pseudo-Hermitian Hamiltonian and 
its Hermitian counterpart, as well as a candidate in experimental realization to simulate the VT 
symmetry breaking. We also show that this model allows conditional perfect state transfer within 
the unbroken VT symmetry region, but not arbitrary. This is due to the fact that the evolution 
operator at certain period is equivalent to the VT operator for the real- valued wavefunction in the 
elaborate PT-symmetric Hilbert space. 
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I. INTRODUCTION 

The transmission of quantum state through a solid 
state data bus with minimal spatial and dynamical con- 
trol is an experimental challenging and a theoretically 
necessary task for implementing a scalable quantum com- 
putation based on realistic silicon devices. S. Bose |l| 
first demonstrated the possibility that in a solid state 
based quantum computer, local interactions can apply 
entangling gates between distant qubits. In principle, 
perfect transfer of quantum state can be implemented 
by specifically engineering chain It is also showed 

that a quantum system possessing a commensurate struc- 
ture of energy spectrum matched with the corresponding 
symmetry can ensure the perfect quantum state transfer 

The aim of this paper is to extend these findings to the 
non-Hermitian system. It is motivated by the interest in 
complex potentials in both theoretical and experimen- 
tal aspects. Much effort has been devoted to establish a 
parity-time {VT) symmetric quantum theory as a com- 
plex extension of the conventional quantum mechanics 
0-Q since the seminal discovery by Bender Q. It is 
found that the non-Hermitian Hamiltonian with simulta- 
neous VT symmetry can have an entirely real quantum 
mechanical energy spectrum and has profound theoret- 
ical and methodological implications. Researches and 
findings relevant to the spectra of the PT-symmetric 
systems were presented, such as exceptional points [l5| . 
s_pectral singularities for complex scattering potentials 
161] , as well as complex crystal and other specific models 
17]. Furthermore, quantum dynamics in open systems 
and PT-symmetric non-Hermitian systems were investi- 
gated 1^- 21 1 . At the same time the VT symmetry is also 
of great relevance to the technological applications based 
on the fact that the imaginary potential could be real- 
ized by complex index in optics [18l. |22h25| ]. In fact, such 
VT optical potentials can be realized through a judicious 
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inclusion of index guiding and gain/loss regions, and the 
most interesting aspects associated with 'PT-symmetric 
system were observed during the dynamic evolution pro- 
cess 

According to the complex quantum mechanics, a non- 
Hermitian Hamiltonian with entire real spectrum can be 
transformed into a Hermitian Hamiltonian via the well- 
established metric operator formalism [30(. Nevertheless, 
the example is rare in which explicit derivations can be 
performed to evaluate the equivalent Hermitian counter- 
part. 

In this paper, we systematically study the VT (parity- 
time reversal)-symmetric non-Hermitian version of a 
quantum network proposed in the work of Christandl 
et al. (HQ. The main results are as follows: (i) The 
non-Hermitian version of a quantum network is exact 
solved and the critical behavior is analytically stud- 
ied at the exceptional point, (ii) We provide an ex- 
plicit mapping between the proposed pseudo-Hermitian 
Hamiltonian and its equivalent Hermitian Hamiltonian. 
The present model is a nice paradigm to demonstrate 
the relationship between the non-Hermitian Hamiltonian 
and its Hermitian counterpart, since all of Hamiltoni- 
ans with different parameters, no matter Hermitian or 
non-Hermitian, are local and have the equally spaced 
spectra, (iii) Furthermore, we investigate the connection 
between the present model and the PT-symmetric hy- 
percube graph, (iv) We also show that such an extended 
model still allows conditional perfect state transfer within 
the unbroken VT symmetry region, but not arbitrary. 
This is due to the fact that the evolution operator at 
certain period is equivalent to the VT operator for the 
real-valued wavefunction in the elaborate PT-symmetric 
Hilbert space. 

This paper is organized as follows. In Section UH we 
present the model and the solutions. In Section IIII1 
we explore the basic properties of the model in its non- 
Hermitian version. In Scction llVl we investigate the met- 
ric and Hermitian counterpart for the extended model. 
In Section |Vj we investigate the connection between the 
present model and the hypercube graph. Section IVII is 
devoted to the dynamics of the model. Finally, we give a 
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summary and discussion in Section [VIII 



II. THE MODEL 

We start with the Hamiltonian for an TV-site tight- 
binding chain with linear potentials, 



iV-l 



H = \ ]T VHN-l) {a\a l+1 + H.c 



1=1 

N r 



(1) 



1=1 



(N+l)-l 



where a\ (ai) is the creation (annihilation) operator at 
site I. For the sake of simplicity we take the units of 
coupling constant as 1. The conclusion of this paper is 
valid for both fermion and boson systems. In the case of 
7 = 0, it is reduced to the model in Ref. 0, H|, which 
has been shown [6| to guarantee that the Hamiltonian 
H(j — 0) evolves states \4>) into V \cj>) at the time 7T, no 
matter what these states are. Here parity operator V is 



given by VajV 1 



ajy +1 _j. This due to the fact that 
the Hamiltonian H (7 = 0) is 'P-symmetric and possesses 
an equally spaced spectrum Q. In the case of nonzero 7, 
the V symmetry is broken. In the following, we will show 
that it still has an equal-spaced spectrum. Here the idea 
is to treat the Hamiltonian as an angular momentum in 
an external magnetic field. 
Defining the operators 



N-l 



j+ = (r) , = Ev / i(M«f 



ai+i, 



(2) 



1=1 



- J+ + J ~ _J + - J 
x r» ' y — 



2i ' 



JV 



Jz = / A — ^ l I ajai 



1=1 



which satisfy the following angular momentum commu- 
tation relations 



=2j z , [r,^] =±j ± . 



(3) 



Then J x ,y,z acts as the angular momentum operator, and 
the Hamiltonian can be rewritten as 



H = J x + 7 J 2 = - ~§ 



where 



t ={B X ,B V ,B Z ) = (1,0,7). 
Obviously, it can be diagonalized as 



(4) 



(5) 



(6) 



with J' n being a n-component of angular momentum op- 
erator, where n = is the unit vector in the field 

direction. In this paper, we concentrate on the single- 
particle invariant subspace, which corresponds to the an- 
gular momentum system with J = (JV — l)/2. Then the 
energy levels are still equally spaced. In this subspace, the 
eigenvector of H can be obtained from that of if (7 = 0) 
by the rotation operator, i.e., 

|V>„) = e -^(7)-/2]./ B (7 = 0)) (7) 
(n = l,2,...JV). 
Here /? (7) is an angle 

/3(7) = arctan (~) (8) 

and \ip n (7 = 0)) is the eigenvector of If (7 = 0), i.e., 

H{l = 0) (7 = 0)) = (^p- - n) |Vn (7 = 0)) 

(9) 
(10) 

(11) 
(12) 



which can be further expressed as 

( 7 = 0))=^d„, (J) |0> . 



Then we have 



where d n ^i (/?) is Winger d-functions 

*m oj) = <r~ l)/2] 03) 

= [(JV - n)\ (n - 1)! (JV - l)\ (I - 1)!] 1/2 

/ „\ N-l+l-n-2v / oX 

(-1)" cosf -sin § 



n-l+2u 



(I- 1 -v)\(N -n-v)\(v + n-l)\v\ ' 

We define ip n (I) = d n .i (/?) to express the orthonormal 
relation as 



^ ipm (I) ipn (0 = &n 



(13) 



It is important to note that the above relation is still true 
for imaginary 7 except the points 7 = ±i. 



III. NON-HERMITIAN PT-SYMMETRIC 
HAMILTONIAN 

Now we consider the Hamiltonian of Eq. ((]]) with imag- 
inary linear potentials by taking 7 — > 17, which can be 
written as 

1 N-l 

H=-^ y/l (JV - (a\a l+1 + H.c.) (14) 



i=i 



z=i 



(JV + 1J-J 
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Mathematically, all the solutions of H can be extended 
to that of % by simply taking 7 — > ij except for 
the points I7I = 1, since it does not induce any sin- 
gularity in the rotation operator e~ lli ^ n ^ Jy . This fact 
accords with complex quantum mechanics. Note that 
the non-Hcrmitian Hamiltonian H is "PT-symmetric, i.e., 
[PT, H] = 0. The antilinear time-reversal operator is de- 
fined as TiT~ l = — i. The phase diagram of this system 
is determined by critical (exceptional) points | -y | = 1. 
This model exhibits two phases: an unbroken symmetry 
phase with a purely real energy spectrum 

£ n = %/T 3 7 5 (^ : -n),n = l ! 2,... ) iV (15) 

when the potentials are in the region [7] < 1 and a 
spontaneously broken symmetry phase with an imaginary 
spectrum when the potentials are in the region I7I > 1. 
In this paper, we only focus on the 'PT-symmetric re- 
gion. Although the spectrum is real, the corresponding 
eigenfunctions are no longer orthonormal with respect to 
the Dirac inner product. One can establish the complete 
biorthogonal set by the eigenfunctions of the Hamiltonian 
TV . Denoting 

\M = E h (0 a \ i°) = E < (0 4 1°) > ( 16 ) 
1 1 

we have 

H f \(j> n ) = e n \<f> n ) . (17) 

Then Eq. (fT3|) leads to the orthonormal relation 

(0„^»)=EC(O<MO = <W (is) 

I 

Furthermore, the complete set obeys the relation 

PT|^) = (-l) n |V„>. (19) 

In the case of 7 = 0, the above relation is reduced to 

V |V« (7 = 0)) = (-1)" |V„ (7 = 0)) . (20) 

It is worth to note that Eqs. (fT9|) and (|20|) have different 
implications: If one obtains a set of eigenfunctions sat- 
isfying Eq. ([9]), they will obey Eq. (l20l) spontaneously. 
However, Eq. (|19p is not necessary for the eigenfunc- 
tions of the Hamiltonian 7^(7). This is due to the fact 
that operator T is a antilinear operator. This issue will 
be elaborated in the next section when we investigate 
the application of Eq. (TT9|) associated with the dynamic 
process. 

Now we consider the case of I7I = 1, which are singular 
points for i\) n (I). However, one can construct the eigen- 
functions of the Hamiltonians W± = % (7 = ±1) as the 
form 

b±>=5>± (o«lio) 

= E(^) 1 ^v / ^( ±l )"" ia 'i )' ( 21 ) 
1 
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FIG. 1. (Color online) Plots of the dressed eigenfuctions as 
functions of for TV = 10 chain. The solid (empty) circle 
indicates the real (imaginary) part of the wavefunction for 
(a) 7 = 0.2, (b) 7 = 0.7, (c) 7 = 0.9, and (d) 7 = 0.998, 
respectively. It can be observed that as 7 tends to 1, all the 
wavefunctions approach to the same function Eq. (|21|l for 
TV = 10, which is plotted in (d) (triangles). Note that the 
scales of the subfigures for various 7 are different. 

where 

V± (l)= (V?) 1 '" JC^^T' 1 (22) 

Acting the Hamiltonians H± on the states \ip±), straight- 
forward algebra shows that 

n± w±) - 0. (23) 

Eigenstates \<p±) are zero-norm states, i.e., 

5>±(0 = 0. (24) 
1 

We can see that I7I = 1 is the boundary of two phases, 
possessing the characteristics of exceptional point: The 
spectrum exhibits square-root-type level repulsion in the 
vicinity of I7I = 1. On the other hand, states \ip±) should 
be the coalescence of the eigenfunctions {|^n)} as I7I — > 
1. To demonstrate this point, we introduce the dressed 
eigenfunctions {IV'n)} by multiplying factor g n (7) to the 
original eigenfunctions {|^n)}, 

<M0 = <M(7)^»(0 (25) 

where 

9n (7) = (V2 cos 1 N (C%-_\y h • (26) 
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We notice that 

lim i, n (I) = <p± (I) , (27) 

±7— vl 

i.e., all the dressed eigenstates coalesce with \<p±) at the 
critical point. For illustration, numerical simulation of 
finite size chains with various values of 7 is performed. 
In Fig. [1] we plot the dressed eigenfuctions including 
real and imaginary parts for N = 10 chain, respectively. 
It shows that all the eigenfunctions tend to \<p±) as 7 
approaches the critical point. 

Finally, we would like to point out that the nature 
of the VT symmetry breaking in such a model is ex- 
clusive. The entire spectrum becomes imaginary and the 
VT symmetry of all the eigenstates is broken at the point 
I7I = 1 simultaneously. The degrees of VT symmetry 
breaking is defined as the fraction of eigenvalues that be- 
come complex p3lj |. the degree of our system is 1. This 
feature should lead to clear signatures in the dynamics 
of the wavepacket. Such a model is a good candidate to 
simulate the critical behavior in experiments. 

IV. METRIC AND HERMITIAN 
COUNTERPART 

Another theoretical interest in the non-Hermitian VT- 
symmetric system is the physical meaning of the Hamil- 
tonian H. When speaking of the physical significance 
of a non-Hermitian Hamiltonian, one of the ways is to 
seek its Hermitian counterparts [32T-|34| possessing the 
same real spectrum. According to the complex quantum 
mechanics, a non-Hermitian 'PT-symmetric Hamiltonian 
can be transformed into a Hermitian Hamiltonian. This 
is achieved by introducing a metric, a bounded positive- 
definite Hermitian operator 77, which can be constructed 
via the eigenstates of V) [30]. However, the obtained 
equivalent Hermitian Hamiltonian is usually quite com- 
plicated [13, H|[ . It is tough to provide an explicit map- 
ping of a pseudo- Hermitian Hamiltonian to its equivalent 
Hermitian Hamiltonian. Fortunately, the present model 
is a nice paradigm to demonstrate the relationship be- 
tween the pseudo-Hermitian Hamiltonian and its Hermi- 
tian counterpart. In this section, we will illustrate the 
basic ideas via the above analytically solution. 

Consider the single-particle case, the metric operator 
X] can be constructed via the eigenstates of H^ as 

n 

= J2 d n,i> G8)aJ|0) (OK 

n,l,l f 

= e -i(p*-p)j y> 

which guarantees the relation 

i-jHT]- 1 = H i . (29) 



It is noticed that the matrix representation of ijH and r\ 
based on the orthonormal basis, say {aj |0)} under the 
Dirac inner product, are Hermitian matrices. Operator 
r\ is called the metric operator since it can be used to 
define the biorthogonal inner product, under which the 
unitary evolution can be obtained. Furthermore, let 

p = s /rj = e-W- p)J * (30) 

be the unique positive-definite square root of rj. Then 
the Hermitian operator p acts as a similarity transforma- 
tion to map the non-Hermitian Hamiltonian % onto its 
equivalent Hermitian counterpart h by 

h = pHp- 1 (31) 

= V0^)H(j = 0). 

It can be regarded as the Hermitian counterpart of the 
non-Hermitian Hamiltonian J? (7). It is fortunate that 
both the non-Hermitian and the Hermitian Hamiltonian 
have simple structure: they possess the localized cou- 
plings and have equally spaced spectra. The physics 
of both the original non-Hermitian Hamiltonian and its 
equivalent Hermitian counterpart are clear: they can 
be regarded as either the lattice model with nearest- 
neighbor couplings (engineered chain with imaginary lin- 
ear on-site potentials) or the angular momentum coupled 
to the external complex field. Additionally, the physics 
of p is also obviously, which presents the operation of two 
successive rotations, e~ 1 ^ / 2 ) J a and e l ^l 2 ^ J y , which de- 
note the rotations about y axis with the angle j3*/2 and 
— /3/2, respectively. 

V. PT-SYMMETRIC HYPERCUBE 

It has been pointed that the one-dimensional system 
H (7 = 0) can be achieved by the projection of the hyper- 
cube Q . In this section, we will extend this approach to 
the non-Hermitian regime and investigate the solutions 
for all the possible projections of the hypercube. 

We start the investigation with an ensemble of non- 
interacting spins in the complex magnetic field if in Eq. 
([5]). The Hamiltonian reads 

d 

-ffspin = ^] Hi (32) 
1=1 

1=1 

where s| is spin-i operator for the lih particle of the d- 
particle ensemble. In z-component spin basis, the matrix 
representation of H spin has the form 

Afspin = ©liMz, (33) 
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FIG. 2. (Color online) Schematics for the system with d = 4 to illustrate the connection among the systems of Eqs. (|14|l . 
(|32[) . and the PT-symmetric hypercube lattice, (a) Four non-interacting spins in the external complex magnetic field, (b) 
The five-column hypercube graph with imaginary potentials, (c) The projection of the hypercube lattice onto the invariant 
subspaces denoted by S. In each invariant subspace, the reduced Hamiltonians are described by Eq. (|14|) with different length, 
respectively. 



where 



VI. PERFECT STATE TRANSFER 



M, 



V) 1 
1 —27 



(34) 



Here © denote Kronecker sum which is defined using the 
Kronecker product ® and normal matrix addition as Mi® 
M £ = Mi ® I2 + I2 ® Mf. . It is different from the direct 
sum of two matrices. Equivalently, such an operation also 
represents the angular momentum coupling of separate 
angular momenta s| and s^. 



According to the graph theory [36(, we notice 
that Mi (7 = 0) also represents the adjacency ma- 
trix of a two- vertex complete graphs fc 2 . Meanwhile 
matrix A/ sp in (7 = 0) is the adjacency matrix of the graph 
Qd = which is constructed by Cartesian product 

of n two- vertex complete graphs k%- Here G1LIG2 de- 
notes the Cartesian product of two graphs Gi and G2, 
and graph Qd is a hypercube graph of d dimensions. 

On the other hand, the total spin S = Y)f—i ~8$ is con- 
servative for the Hamiltonian -ff S pin, i-e., [S 2 , i? S pm] = 0. 
Then matrix M sp i n can be diagonalized in each invariant 
subspace denoted by S. In each invariant subspace with 
S = d/2, d/2 — 1, (or 1/2), the matrix is the same as 
that of Eq. (fT4|) with N = 2S + 1 in single-particle case. 
For the case of N = d + 1, it accords with the analysis 
for the engineered chain in Ref. Q. The column in the 
hypercube graph corresponds to the set of the states with 
the same z-component spin S z = 53/=i s i ■ I n Fig- H]wc 
take d — 4 for an example to illustrate our analysis. 



Now we turn to the dynamics of the non-Hermitian 
system. We notice that the imaginary linear potentials 
break the V symmetry but retains the VT symmetry. 

In the following, we will explain how to use such system 
to realize the perfect state transfer. The time evolution 
of the initial state |$ (0)) can be expressed as 



|*(t)) = e- 

where C n = (cf> 
(TT51) . At time t 
(TT9l) we have 



-mt 



l$(o)) = ^a 



-iE n t 



\lpT, 



(35) 



|$(0)) can be obtained from the Eq. 
r = 7r/-\/(l — 7 2 ), from Eqs. ([IB]) and 



|$(r))=^C„(-in^) 

n 
n 

n 

For the state with all real C n , we have 

I$(t)) = :pt !<£(())}, 



(36) 



(37) 



i.e., at the instant r, the time evolution operator acts 
as VT operator. For 7 = case, Eq. (|3"T|) reduces to 
l$( r )) = for arbitrary states |$(0)). As we 

mentioned in the last section, Eq. (IT5)) is not necessary 
for a given eigenstate. Nevertheless, Eq. (|37|) can pro- 
vide a way to construct the state to satisfy the Eq. (|19l) . 
In the following we will exemplify this point and its ap- 
plication. 
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FIG. 3. (Color online) Plots of the time evolution for real-valued quantum states in N = 100 chain. The initial state has the 
form of Eq. <|41[) . The black (red) line indicates the real (imaginary) part of the wavefunction for (a) 7 = 0.05, / = 1, (b) 
7 = 0.05, I = 25, (c) 7 = 0.1, / = 1, (d) 7 = 0.05, I = 100, respectively. It shows that the width of the initial state becomes 
wider as I or 7 increases. We can see that the final state at the time t = r is the VT function of the initial state. In (a, b, c) 
the amplitudes of the evolved state increase during the internal [0, r/2], then decrease during [r/2, r]. In contrast, it decreases 
during [0, r/2], then increases during [r/2, r] in (d). Note that the scales of the subfigures are different. 



With the same mechanism of the action of operator T 3 , which can be expanded with eigenstates 
operator VT also takes the role of perfect quantum state 
transfer. The flaw of this scheme is that it only applicable 

for some specific state. However, if there exists local ft t |g^ _ \ip m ) (cf> m \a? |0) (38) 

state satisfying Eq. (|3T|) . it has a potential for future ^ 
applications. r w-1 1 

= £4? J (/3)l^n). 



In the following, we will show that the local state can 
be constructed to perform perfect state transfer for small 
|7|. Consider a local state at the end of the chain a\ |0), 



Obviously, such a state does not satisfy Eq. (|3"T)l and 
cannot be perfectly evolved to the state ajy |0). Other- 
wise, one can construct a state satisfying Eq. (|3"T|) based 
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on the state a[ |0) in the way 

where Qi = 2 + (l — 7 2 ) 1 W is the normalization factor. 
Note that the expansion coefficients of |1) are all real, so 
that it can evolve to the state VT\1). Now we will prove 
that such a state is local in the case of I7I <C 1. Actually, 
rewriting the state 1 1) in the basis {aj |0)} and neglecting 
the high-order terms of the Taylor expansion, we obtain 

|T> =a t 1 |0)-z^ v ^V^T)4|0), (40) 

which is a local state. At time t = r = 7r/y(T— 7 2 ), it 

evolves to state ^ajy + ijy/ (N — 1) /2a] v _ 1 ^ |0). 

In order to demonstrate and verify the above analysis, 
we perform the numerical simulation for a finite TV-site 
system. The initial wave function has the form 

10 = Tff £ { d Cf 1 + d iT ] (P)} Mm) , (41) 

r w-1 1 r w-i i 
where Q t = 2 + 2 £^ =1 4 j J (/?) 4,; 2 J (0) is the 
malization factor. The evolved wave function 

$(n,t) = {0\a n e- im \T) (42) 

is computed via exact diagonalization method. We plot 
the real and imaginary parts of $ (n, i) as function of 
time in Fig. [31 respectively. It shows that the evolution 



process is different from that in a Hermitian system: the 
Dirac inner product of the evolved state is not conserva- 
tive. The time evolution during intervals [0, r] and [r, 2t] 
are completely different processes. 



VII. SUMMARY AND DISCUSSION 

In this paper, we have shown that adding the VT- 
symmetric potentials on the well-studied Hermitian 
quantum network 0, 0] constructs a exactly solvable non- 
Hermitian model which allows conditional perfect state 
transfer within the unbroken VT symmetry region, but 
not arbitrary. This model has applicability and rele- 
vant to the physical situations since there exist local 
states which can be transferred perfectly across long dis- 
tance. In the theoretical aspect, this work provide a 
nice paradigm to demonstrate the relationship between a 
pseudo-Hermitian Hamiltonian and its Hermitian coun- 
terpart in the framework of the complex quantum me- 
chanics. On the other hand, the simultaneity of the onset 
of the VT symmetry breaking for the whole eigenstates 
should lead to remarkable phenomena in the dynamics 
of a wavepacket. This result suggests the evident obser- 
vation of the VT symmetry breaking in optical system 
with complex index. 
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